Calibration of optical tweezers with positional detection in the back-focal-plane 
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We explain and demonstrate a new method of force- and position-calibration for optical tweez- 
ers with back-focal-plane photo detection. The method combines power spectral measurements of 
thermal motion and the response to a sinusoidal motion of a translation stage. It consequently 
does not use the drag coefficient of the trapped object as an input. Thus, neither the viscosity, nor 
the size of the trapped object, nor its distance to nearby surfaces need to be known. The method 
requires only a low level of instrumentation and can be applied in situ in all spatial dimensions. It 
is both accurate and precise: true values are returned, with small error-bars. We tested this exper- 
imentally, near and far from surfaces. Both position- and force-calibration were accurate to within 
3%. To calibrate, we moved the sample with a piezo-electric translation stage, but the laser beam 
could be moved instead, e.g. by acousto-optic deflectors. Near surfaces, this precision requires an 
improved formula for the hydrodynamical interaction between an infinite plane and a micro-sphere 
in non-constant motion parallel to it. We give such a formula. 
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I. INTRODUCTION 



In order to use optical tweezers as a quantitative in- 
strument for position and force measurements, the detec- 
tion system must be calibrated. One calibration method 
interprets the power spectrum of thermal Brownian mo- 
tion of a trapped object P, S Another calibration 
method interprets the displacement of a trapped object 
in response to a known flow past it 0, Q ■ Both these 
methods require that the drag coefficient of the trapped 
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object is known. Here we combine the two methods into 
one. Then there is no need to know the drag coefficient. 

Mathematically speaking, the combined method mea- 
sures an extra quantity that allows the elimination of 
the drag coefficient from the calibration procedure. This 
has some advantage: Methods that use the drag coef- 
ficient as input, use Stoke's law to calculate it, hence 
rely on assumptions about the object's shape and radius, 
and about the viscosity of the surrounding fluid. This 
contributes to the error on the final calibration, as do 
hydrodynamical interactions with nearby surfaces jfl]. 

With the combined method presented here, these 
sources of error have been eliminated. The drag- 
coefficient of the trapped object can in fact be mea- 
sured with the present method, and this with precision, 
as shown below. Finally, the method presented here is 
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simpler to implement and requires less instrumentation 
than other methods that use additional lasers 0, or 
AODs 



II. MATERIALS AND METHODS 

Measurements were done with two different optical 
tweezers systems. Both use the trapping laser for po- 
sition detection in the back focal plane. One system has 
a long working distance and was used to test the method 
far from surfaces (30 /im): This experimental setup is 
described in detail in and briefly here. The instru- 
ment is based on a custom-built inverted microscope with 
a Nikon, 60x, 1.2 NA, 0.2mm working distance, Plan- 
Apo, water immersion objective. The laser is a 1064 nm, 
Nd:YAG (Spectra-Physics Millennia IR). Position detec- 
tion is done with a position-sensitive photodiode (UDT 
DLS-20). Flow-cells with a volume of 10 /il (dimensions 
8 X 20 X 0.06 mm'^) were assembled by placing a covcr- 
slip on top of a microscope sHde separated by spacers of 
double-sided sticky tape. A dilute solution of beads was 
flowed in, and the ends were sealed with nail polish to 
avoid sample evaporation. The flow-cell was mounted up- 
side down on a Physik Instrumente piezo-electric trans- 
lation stage (P-527.2 CI). 

The other system has a short working distance and 
was used to test the method close to surfaces (0-3 um): 
This experimental setup is described in detail in [ll|. 
Briefly, it consists of a modified Zeiss Axi overt 135 TV 
microscope equipped with a Zeiss, lOOx, 1.3 NA, Plan- 
Neofluar, oil-immersion objective. The laser is a 1064 nm, 
Nd:YV04 (Smart Laser Systems GmbH, Berlin, Ger- 
many). Position detection is obtained with a standard 
quadrant photodiode, QP50-6SD (Pacific Silicon Sensors 
Inc). Signals were recorded with a 24-bit data acquisi- 
tion card (NI 4472, National Instruments) which has a 
45 kHz alias-free bandwidth. Temperature was measured 
with type-T thermocouples, accurate to within 0.1 °C 
(IT-23 Physitemp, Clifton, NJ, USA). Flow-cells with a 
3 mm- wide channel were assembled by placing a 18mm^ 
cover-slip on top of a No. 1.5 22 mm^ cover-slip sepa- 
rated by a layer of parafilm. The parafilm was melted by 
placing the sample on a 100°C hot-plate. Cooling then 
glued the cover-slips together. A dilute solution of mi- 
crospheres was flowed in, and the ends were sealed with 
vacuum grease to avoid sample evaporation. The flow- 
cell was mounted on a Physik Instrumente piezoelectric 
translation stage (P-733.2 CL) with built-in capacitance 
position detection pre-calibrated to within 0.1%. For 
measurements close to the surface, the cover-slips were 
treated as described in to reduce the influence of 
surface potentials. 

The following beads were used: Silica beads from 
Bangs Laboratories (9025 Technology Drive, Fishers, IN 
46038-2886, USA), catalog code SS04N, lot number 5303, 
were 1.54 /xm in diameter with a coefficient of variation 
(standard deviation divided by the mean value) of 10% 



listed by the manufacturer. Polystyrene microspheres 
from Polysciences (Warrington, PA 18976, USA), cat- 
alogue number 07307, lot number 50602, were 528 nm 
in diameter with a 2% coefficient of variation listed the 
manufacturer. Transmission electron microscopy (TEM) 
of the latter microspheres indicated a 1.2% coefficient of 
variation. 



III. THEORY 

In the following we present the theory for a calibration 
procedure in which the flow-cell is moved sinusoidally a 
known distance relative to the trapping laser by a piezo- 
electric translation stage. If, instead, the trap is moved 
a known distance relative to the flow-cell, the same for- 
mulae apply, provided the detector is placed in the back 
focal plane of the condenser. 

In order to keep the presentation as simple as possi- 
ble, we use Einstein's simple theory for Brownian motion. 
The presentation given here carries over unchanged to the 
full, hydrodynamically correct theorVjincluding possible 
filters, electronic and/or parasitic |2,Ql- We present the 
relevant formulas in Appendix IHl 

A. Equation of motion 

A microsphere suspended in water is trapped with op- 
tical tweezers inside a flow-cell. The stage moves the 
flow-cell sinusoidally relative to the optical trap with a 
frequency /drive and an amplitude A, see Fig. ^ while 
the trap remains at rest in the laboratory system. The 
position of the stage as a function of time t is 

a;drivo(i) = ^sin(27r/drivci) ■ (1) 

The stage velocity Udrivc(0 = idrivc(i) also corresponds 
to the velocity of the water in the flow-cell far away from 
the bead, since the water is at rest relative to the flow-cell 
(see Appendix IbI). Ignoring hydrodynamical and incrtial 
effects, the Langevin equation of motion for a spherical 
bead in the trap is 

^drivc 

(t)] + Kx{t) = Frit) (2) 

where x{t) is the position of the bead relative to the cen- 
ter of the trap, 7 is the drag coefficient, and k is the 
trap stiffness. The first term on the left-hand-side is the 
drag force, which is proportional to the velocity of the 
bead relative to that of the stage. The second term on 
the left-hand-side is the trapping force. The right-hand 
side is the random thermal force driving the Brownian 
motion. It is assumed to have the statistical properties 
of white noise. 

Frit) = = -/V2D({t) (3) 

where k^T is the Boltzmann energy at absolute temper- 
ature T, the diffusion coefficient is given by Einstein's 
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FIG. 1: Positions of the piezo-stage (a) and the trapped 
bead (b) for a stage moving sinusoidally with frequency 
/drive = 32 Hz and amplitude A — 150 nm in the x-direction. 
(a): Left, time series of stage position. Right, histogram of x- 
coordinate of the stage position. The sinusoidal movement re- 
sults in two clear maxima, (b): Left, time series of the bead's 
a;-coordinate in volts, as given by the signal from the pho- 
todiode. The amplitude of the sinusoidal response is smaller 
than the amplitude of the thermal motion. Consequently the 
sinusoidal shape is masked by the Brownian motion of the 
bead, and the maxima associated with the sinusoid disappear 
in the histogram of visited positions (right). 



relation D = k-sT/^, and ^ is the normalized white noise 
expressed with Dirac's delta function S, 



(4) 



Here (...) denotes the expectation value that results from 
averaging over the thermal noise. 

B. Solution to the equation of motion 

Since the equation of motion Q is linear with two force 
terms. Ft and 7Wdrivc7 its general solution can be written 
as a sum of two terms, one for each force, 



X{t) = XT{t) + a;rosponsc(i) , 

after transient initial behavior has died out. Here 



(5) 



xrit) = V2l5 f dt'e-2-/c(*-t')^(i') ^ (6) 

^ — oo 



"^response 



(t) = 



+ (/c//drivc)2 



(7) 



where we have introduced the corner frequency fc = 
k/{2ttj), and tiag = [arctan(/drivc//c) - 7r/2]/ (27r /drive)- 
Figure ^3 shows an example of an experimentally deter- 
mined trajectory x(t) of a bead in a trap, as described in 



Eq. (ISJ. In this case, the stochastic thermal motion dom- 
inates and almost hides the driven, deterministic compo- 
nent of the motion. 

In principle, we now could calibrate by fitting 
2^rosponso(^) in Eq- to data like those shown in Fig.^. 
However, this is not a reliable procedure Instead, 
we Fourier transform theory and data to the frequency 
domain where parameters are determined with optimal 
precision because the theory is simpler there. 



From Eqs. (|5] 

of x{t) is 



C. Power spectrum 



it follows that the Fourier transform 



xif) = / dte"-/*a;(i) 

J — OO 

^ iif)V2D 
Mfc - */) 



(8) 



2*v/l + (/c//drivc)^ 



m + /drive) - Sif - /drive)] 



where ^(/) is the Fourier transform of ^(t). Consequently, 
the expectation value for the one-sided (/ > 0) p ower 
spectral density (PSD) of the bead positions is [13 

2(ii(/)r 



pif) 



PTif)+Pr. 



if) 



(9) 



D 



t ' 2t -C2 < t2\ ' 1 I -f2 / -f2 /drive) 

W-»oo TT^if'^+f^) l + /cV /drive 

where inisr is the measurement time (see also |2, Eq. (8)]). 
This PSD consists of the familiar Lorentzian (first term, 
Pt), plus a delta-function spike (second term, Presponso) 
at the frequency with which the stage is driven. The 
Lorentzian originates from the Brownian motion of the 
bead in the parabolic trapping potential, and is hereafter 
referred to as the "thermal background." 



IV. HOW TO CALIBRATE 

Experimentally, positions are measured in volts, a;™'*. 
Assuming linearity 0, 0| , 



a;(t) =/3a;™"(i) , 



(10) 



positions are known in meters once the calibrations fac- 
tor P has been determined. It can be determined from 
the measured PSD: From Eq. (|10|) follows that the ex- 
perimental PSD p™" is measured in (volts) ^ x (seconds) 
and 



-Prosponse(/) — /3 -Presponse (/) 



(11) 



4 



Here, -Pr™ponsc known experimentally and, as can be 
seen from Eq. so is -Prcsponsej since A and /drive are 
known a priori, and fc is known experimentally. So the 
calibration factor (3 is the only unknown in Eq. 
hence determined by this equation. 

In practice the experimental measurement time is al- 
ways finite. This is taken into account in the expression 
for P given below. 

The procedure described here works in bulk and near a 
surface. If used near a surface, the resulting calibration is 
specific to the distance from the surface. So is the value 
found for the drag coeSicient. Faxen's law (Eq. I)D5|I 
below) and formulas in [TllIT^ Sect. 7.4] relate the drag- 
coefficient to its bulk value for the lateral and axial di- 
rection, respectively. 



A. Positional calibration 

The desired calibration factor for distances is 

P = y/Wtb/W,^ , [/3] = m/V , (12) 

where Wox is the experimentally determined power in the 
spike, measured in (volts)^, and Wth is the same quantity 
measured in (meters)^. The latter follows from Eq. @, 

/•/Nyq 

Wtl. = / Prcsponsc(/)d/ = , ^ % . , (13) 

-'0 ^ T" /c / /drive 

whereas Wcx is found by observing that in an experiment 
the measurement time inisr is not infinite, so the spike is 
no longer a delta-function. As detailed in Appendix lUl 
its height is now proportional to inisr, and its width is also 
finite in general. However, if frnsr is an integer multiple of 
the period of the stage movement, then the spike consists 
of a single datum [l^ , see Fig. |21 and consequently 

W^ex = (P™" (/drive) ~ (/drive)) A/ , (14) 

where -P™'* (/drive) is the height of the spike, i.e., the ex- 
perimentally determined value of the PSD at /drive and 
P™'* (/drive) is the PSD of the thermal background at 
/drive. Both are measured in (volts)^ x (seconds). The 
width A/ = l/trnsr corrcsponds to the frequency resolu- 
tion of the PSD. 

The thermal background in this expression can be 
treated with three different levels of precision: (i) Do 
not subtract it at all, if it is truly negligible compared 
to the height of the spike, (ii) Interpolate its value from 
a plot of the power spectrum. This approach can also 
be taken if the PSD is contaminated with low-frequency 
power, that is external to the model used here, e.g. from 
electronic noise or the pointing instability of the laser 
beam, (iii) Use the value at /drive of the theoretical ex- 
pression for the thermal background P™'*(/) after it has 
been fitted to the experimental thermal background. 




Frequency (Hz) 

FIG. 2: Power spectrum of a 528 nm diameter polystyrene 
bead held in the laser trap with a corner-frequency fc = 
2065 ± 5 Hz. The sample moves sinusoidally with amplitude 
A = 150 nm and frequency /drive = 32 Hz. The power spec- 
trum shown is the average of 100 independent power spec- 
tra. It consists of a thermal background caused by Brownian 
motion, plus a spike at /drive- The sampling frequency was 
/sample = 65 536 Hz, the measurement time for each spectrum 
imsr = 1 s, and the temperature 24.4° C. For clear illustration, 
the measurement time used here is eight times longer than the 
one we typically use for calibration. 

B. Force calibration 

The trap's force on a bead is kx. With the position 
detection system calibrated, the displacement x is mea- 
sured in meters. To determine the trap stiffness k, we 
use the definition of fc and Einstein's relation to write 
K = 2t: fckBT/D. Next, we determine /c and Z)™'*, the 
diffusion constant in (volts)^/(seconds), by fitting the 
first term in Eq. Q to the thermal background in the 
experimental PSD. Using D ~ f3^D^°^^, we arrive at a 
practical formula for force-calibration 



All variables on the right-hand side are measured experi- 
mentally in the calibration process, including Wth, which 
is determined as the right-hand side of Eq. The 
local temperature T of the liquid can typically be deter- 
mined with sufficient accuracy by direct measurement in 
or near the flow-cell: T is the absolute temperature, so 
an absolute error of 1 K results in a relative error of only 
0.3%. 

As a side effect of this calibration procedure, we find 
an experimental result for the drag coefficient, 

where again all variables on the right-hand side are mea- 
sured experimentally in the calibration process. 
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Use of the hydrodynamically correct theory for Brown- 
ian motion changes Eq. H12|) only by replacing Wth with 
hydro (sce Appendix IDjl . The values found for /s^cx 
and 7ox with these equations do change, however, be- 
cause more correct values for fc and D™'* result from 
fitting the thermal background with the hydrodynami- 
cally correct theory. The value found for Wcx is also 
affected a little through the subtraction of the thermal 
background in Eq. (|14|l . 

V. EXPERIMENTAL RESULTS 

In order to test the accuracy and advantage of the cal- 
ibration method described above, we here compare our 
experimentally determined values for the drag coefficient 
with values calculated from vendor information. Deep 
in bulk, 21 silica beads with a diameter of 1.54 /^m were 
studied with an optical tweezers system with a water- 
immersion objective. At various distances close to the 
cover slip (< 3/im), 24 polystyrene beads with a di- 
ameter of 528 nm were studied with an optical tweezers 
system with an oil-immersion objective. Fitting of the 
PSD to the thermal background recorded for each bead 
at each of its positions considered was done with the 
highest possible precision, using either published Mat- 
Lab routines [la, [l3 or custom written software in Lab- 
View ^3 . These fitting routines take into account hydro- 
dynamics corrections, aliasing, parasitic filtering in the 
photo diode, and electronic filters in the data acquisition 
system 

A. Measurements in bulk 

First we measured the drag coefficient far from sur- 
faces, where 7ex can be compared directly to Stokes's for- 
mula 7o = GnriR, with the bead diameter 2R = 1.54 /xm 
taken from the specifications of the producer, and the 
viscosity rj = 0.93mPa-s calculated for water at the mea- 
sured temperature, 23.0 °C. Beads were trapped near the 
bottom of the flow-chamber (silica beads are heavier than 
water) and brought to the middle of the flow-cell, 30 ^m 
from the bottom and the top, to minimize the effect of 
nearby surfaces. The x-axis — the direction of motion — 
was chosen perpendicular to the long axis of the flow- 
cell and perpendicular to the direction of the incoming 
laser light (the z-axis). The experimental parameters 
were A = 208 nm and /drive = 28 Hz. For each of the 
21 beads, fc, D™'*, and Wcx were determined from its 
power spectrum recorded in bulk. 

The average of values measured for 7ex was 7ex = 
( 13.4 ± 0.2) nN-s/m (mean ± standard error on the mean 
(SE), n = 21 beads), and the SD of measured values for 
7cx was 0.8 nN-s/m, i.e., the measured drag coefficient 
showed a coefficient of variation of 6%, which agrees with 
the 10% variation in bead diameter listed by the manu- 
facturer because of our limited sample size. The value 



for the average drag coefficient expected from the man- 
ufacturer's information is (13.5 ± 0.3) nN-s/m (mean ± 
SE) where the SE was obtained by adding the errors as- 
sociated with viscosity and variation in bead diameter in 
quadrature divided by n. These results suggest that our 
calibration method is both accurate and precise. 

B. Measurements near surfaces 

In practice, experiments are often done close to a sur- 
face, such as a cover slip. Working close to a surface — 
in our case < 3 fim — complicates standard calibration 
techniques, especially if a high numerical aperture oil- 
immersion objective is used. In that case aberrations 
arising from a refractive index mismatch at the glass- 
water interface cause a linear decrease in stiffness away 
from the surface and a focal shift 0, 0, |23| ■ 

1. Acquisition and fitting 

The experimental parameters were, A = 150 nm, 
/drive = 32 Hz, /sample = 65 536Hz, t^sr = l/8s, and T = 
24.4°C. The beads used had, according to the producer, 
a diameter of 2R = 528 nm with a 2% coefficient of vari- 
ation. Power spectra that resulted from averaging 100 
independent spectra were fitted with a custom-written 
least-squares fitting routine implementing a Levenberg- 
Marquardt algorithm (Labview, NIV. Each datum was 
weighted by its theoretical error bar H- The propagated 
error on the fit parameters was calculated as the square 
root of the diagonal elements of the covariance matrix 
multiplied by the reduced x^-value. The fit was done in 
the frequency interval [8 Hz : 25 kHz] , omitting the single 
datum at the stage frequency /drive — 32 Hz, and using 
Eq. ljD2|) for the thermal background, with Eq. (35) in 
describing parasitic filtering. 

For every trapped bead the calibration factor f3, the 
trap stiffness Kqx, and the drag coefficient 7ex were found 
at each of 50 distances from the surface. The exact sur- 
face position was obtained from a fit of Faxen's law to 
the measured drag coefficients, taking the focal shift into 
account; see Appendix IdI and .HJ. Thus, the quoted val- 
ues for 7ox are the measured values extrapolated to bulk, 
hence directly comparable to 70. 

2. Temperature 

We measured the temperature with a small thermo- 
couple introduced into the flow cell, while simultaneously 
recording the room temperature and the temperature of 
the imaging and condenser objective. The temperatures 
of these two objectives differed by 0.5 °C. The tempera- 
ture inside the flow cell was intermediate between these 
and was measured to within 0.2 °C. We estimated an up- 
per limit of 0.5 °C for the local temperature increase due 
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to laser heating |2l|. Thus, the propagated error on 7ox 
in Eq. H16|l would be less than 0.2% if uncertainty about 
the temperature were the only source of error. 



Drag coefficient 



than this poly-dispersity. The coefficient of variation of 
the bead-population's various 7ox-values was 1.5%, which 
is comparable to the 2% coefficient of variation for the 
bead diameter listed by the producer, and the 1.2% mea- 
sured independently by us using TEM. 



Figure 13 shows the experimentally determined values 
for the drag coefficient, 7cx, here given in units of the 
theoretically expected drag coefficient, 70, for 24 individ- 
ual beads. The error-bars on the individual data-points 
are the propagated errors from the fit to Faxen's law 
and the uncertainty in the temperature. This high pre- 
cision results from the long measurement time and the 
50 determinations of 7ex for each bead. Any systematic 
error — e.g. an undetected error in the specifications of 
the piezo-stage — will offset the mean value of 7cx, but 
will not change the position of the data-points relative to 
each other. In other words, systematic errors directly in- 
fluence the accuracy of the method, but do not influence 
the precision. We check the accuracy by comparing 7cx to 
7o, including its estimated errors. The shaded area shows 
the 2.3% error on 70 from the propagated uncertainty on 
the viscosity (temperature) and the bead radius. 



1.05 



1.00 



0.95 



§ 5§ 

- ^e- - 



A' 



10 15 
Bead number 



20 



25 



FIG. 3: Results from the calibration of 24 individual beads 
trapped close to a surface. Circles show the measured drag 
coefficients Eq. (Ilfcil , extrapolated to their bulk values accord- 
ing to Faxen's law, in units of the theoretically expected drag 
coefficient in bulk. Error-bars are the propagated errors from 
the fitting routines used. The horizontal solid line and the 
two dashed lines denote the mean (7^) ± SD of the 24 mea- 
surements. The accuracy of the method is illustrated by the 
agreement of the two values 7^/70 = 1.007 ± 0.003 (mean 
± SE). The shaded region denotes the estimated 2.3% uncer- 
tainty on 70. 

The average of the measured values for the 24 beads, 
TW, was 1.007± 0.003 (mean ± SE, n = 24) in units of 70. 
As seen in Fig.|31 the experimental value for an individual 
bead may differ several percent from the theoretically 
expected value 70, even if the average value does not. 
Thus, if 7o is used for calibration, stochastic errors of 
several percent are expected due to the poly-dispersity 
of the bead radii. 

We were able to measure the poly-dispersity of the 
bead population because the precision of our calibration 
method (error-bar on single-bead datum) was smaller 



VI. DISCUSSION 

The main strength of the calibration method presented 
here is that it is independent of estimates regarding the 
drag coefficient and that it can be performed on the loca- 
tion of the experiment. It shares this strength with other 
methods 0, Q , but unlike those it does not require two 
laser beams or acousto-optic deflectors. 

The method presented here has some resemblance to 
a method used in [2^ There, a bead embedded in 
a gel was driven sinusoidally through the detection laser, 
and the resulting spike in the PSD was used to estimate 
the sensitivity of the detection system. However, that 
method was not in situ since the bead used for the cali- 
bration could not be used for later experiments. 

In Appendix we discuss a number of issues associ- 
ated with the use of optical tweezers in general and our 
implementation of the calibration method in particular. 



VII. CONCLUSIONS AND OUTLOOK 
A. Recommended approach 

We suggest the following steps in experiments using a 
piezo-electric translation stage: 

1. Trap the object of interest at the position of inter- 
est. 

2. Drive the stage at any frequency /drive that does 
not excite resonances in the system. We typically 
used 16 Hz or 32 Hz, low enough that the fluid 
moves with the stage, see Appendix IbI 

3. Collect position data from the photo detection sys- 
tem and the piezo-stage for a time t^sr that is an 
integer number of stage periods — i.e., tmsr/drivc = 
integer — to avoid leakage, and with a sampling 
frequency that makes the number of data points, 
^msr/sampie, & powcr of two, SO thc PSDs of posi- 
tions can be obtained by FFT. 

4. Determine A and /drive from thc PSD of thc piezo- 
stage positions (or use previously calibrated val- 
ues). 

5. Determine and /c from a fit to the thermal 
background in the PSD of the bead-positions. 

6. Determine the value of the spike in the PSD at 
/drive and calculate /3 and k using Eqs. (|12|l and 
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In Point 12 the suggested values for /drive are powers 
of two for ease of calculation of the appropriate ty^si and 
/sample- In Point El we recommend using the full theory 
for fitting the PSD if less than 10% systematic error is 
desired. The relevant equations a re g iven in Appendix IdI 
and in 0,01) and software in p^lrnj. 



B. Calibration of the axial direction 

The total intensity of the laser light that impinges on 
the position detector in the back-focal plane, contains in- 
formation about the displacement of the trapped object 
in the axial direction, i.e., in the direction of the optical 
axis, the z-axis, It is generally possible to calibrate 
this dimension by moving the sample relative to the laser 
along the z-axis. If acousto-optic deflectors are used to 
move the laser relative to the trapped object 0], this is 
impossible. However, for spherical objects the bulk drag 
coefficient is the same in all directions, and it is suffi- 
cient to determine the drag coefficient in one dimension 
in order to calibrate all three dimensions. When work- 
ing close to a surface, the drag coefficient depends on the 
direction of motion. However, if the bead is spherical 
and the distance to the surface is known or determined 
as, e.g. in Subsect. I VBL it is again sufficient to determine 
the drag coefficient in just one direction. Faxen's law and 
formulas given in lUlllJ, Sect. 7.4] can be used to correct 
for the distance dependence of the drag coefficient for the 
remaining lateral and/or axial dimensions, respectively, 
depending on which dimension was previously calibrated. 



C. Imaging versus non-imaging positional detection 

Here we demonstrated our method using positional de- 
tection in the non-imaging back-focal-plane of the con- 
denser. In the back-focal plane an interference pattern is 
detected. This pattern arises as the laser light scattered 
from the bead interferes with the unscattered laser light. 

Position detection in the image plane should also work, 
if the laser-trap remains stationary in the laboratory co- 
ordinate system, while the flow cell is driven. For ex- 
ample, an image of the bead could be projected onto a 
photodiode or a camera. With such a setup, calibration 
can also be obtained by simply moving the detector a 
known distance relative to the image while recording the 
response 25^. 

Another approach moves the trapping laser relative to 
the flow-cell, e.g. using acousto-optic deflectors, galvano- 
mirrors, or some other beam-steering apparatus This 
approach works, if back-focal-plane detection is used: A 
pure translation of the laser in the image-plane produces 
no signal in the back-focal-plane, only a motion of the 
trapped object relative to the laser is detected. An ad- 
vantage of this approach is that open samples can be 
used, because the sample is not moving. 



D. In situ measurements 

The method presented here is implemented strictly in 
situ. Therefore, it should be applicable in situations 
that have so far eluded accurate measurements of posi- 
tions and forces, e.g. when trapping spherical structures 
of unknown refractive index and size in the interior of 
cells. The method described here could also be used in 
microfluidic lab-on-a-chip devices, for measuring forces, 
viscosities, or temperatures inside micron-sized channels. 
Generally, position and force calibration should be pos- 
sible in an arbitrary geometry, because we do not need 
to know the corrections to the drag-coefficient due to the 
proximity of surfaces in the approximation used in Eq. (0) 
where the frequency-dependence of the drag-coefficient 
is neglected. More generally, all that is needed for the 
present method to apply with precision is a system that 
displays a linear response to forces, as in Eq. ©I The 
drag should be proportional to the trapped object's ve- 
locity, but may depend on frequency as in Appendix ^ 
as long as the functional form of the frequency depen- 
dence is known, or can be modelled, and the trapping 
force should be Hookean. 



VIII. ACKNOWLEDGMENTS 

HF and FJ thank the Isaac Newton Institute for Math- 
ematical Sciences and its program Statistical Mechanics 
of Molecular and Cellular Biological Systems for hospi- 
tality while this work was initiated. HF thanks Steve 
Block for the opportunity to present an early version of 
the present work at his Winter Workshop on Biophysics, 
The Biophysics of Single Molecules, Aspen, CO, January 
2-8, 2005. SFT-N thanks Marco Capitanio for technical 
assistance. We all thank the anonymous reviewer for a 
thoroughly critical effort that improved the presentation 
a good deal. 

APPENDIX A: EXPERIMENTAL ISSUES 

We considered the following issues in the implementa- 
tion of our calibration method. 

1. Stage response. Piezo stages have a finite response 
time, but this is not a problem here, since we drive 
with a sinusoidal signal. This is the signal, which 
another periodic signal will degenerate towards, if 
its period is shorter than the response time of the 
stage. Also, any deviation from the chosen signal 
shows up in the experimental power spectrum as 
distinct, isolated higher modes at frequencies that 
are integer multiples of the driving frequency. 

Mechanical resonances of the experimental setup 
can be exited. Our setups have resonance frequen- 
cies starting at approximately 400 Hz. The excita- 
tion of a resonance will give rise to additional power 
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in the PSD at the resonance frequency, hence is 
under experimental control. Mechanical crosstalk 
between axes occurs, as seen in Fig. ^ However, 
as the equations of motion are linear, this does not 
influence the calibration. 

2. Photodiode response. Photodiodes may act as fil- 
ters, but the effect is well understood and can be 
measured and accounted for "2^. Depend- 
ing on the type of diode used (quadrant or position 
sensitive), the linearity of the response to bead dis- 
placements may also vary. The region of linearity 
is easily found by moving a stuck bead through 
the laser focus. The photodiode must be aligned 
with the piezo-stage. This is achieved by driving 
the stage along one of its axes while rotating the 
diode until the power at the driving frequency is 
maximized for the corresponding axis. 

3. Crosstalk. Crosstalk between the axial (z) and 
the lateral {x and y) channels from a quadrant- 
photodiode may lead to underestimates of /c of 
up to 10%. This crosstalk shows up as additional 
power in the PSD of the x, y-channels. Since the ax- 
ial corner frequency is smaller than the lateral ones 
pyf. crosstalk raises the plateau of the lateral sig- 
nals. A likely source of the crosstalk is differences 
in amplification of the signals from the diode's four 
quadrants. By repositioning the diode relative to 
the laser so that a small offset from the center po- 
sition (i.e. zero volts in x and y) is introduced, the 
crosstalk can be minimized by maximizing the lat- 
eral corner frequencies that are returned from fits 
to the PSD '3,flH. 

4. Laser heating of the liquid. The laser heats up the 
liquid locally, resulting in a decrease in viscosity. 
We looked for this effect by varying the laser inten- 
sity, but did not find any such effect when working 
close to surfaces. This result is to be expected be- 
cause the glass cover-slip acts as a heat-sink |2ll |. 

5. Hydrodynamic response of the sample. The calibra- 
tion method presented here is conditioned on the 
liquid co-moving with the stage. In Appendix IBl we 
calculate the response of the liquid to the oscilla- 
tory movement of the stage. Close to the surface, 
the no-slip boundary condition entrains the liquid. 
Further into the sample, the degree of entrainment 
depends on the height of the sample d and the drive 
frequency. If d^TT /drive < ~ Imm^/s, where v is 
the kinematic viscosity, the liquid co-moves in the 
entire flow-cell. 

6. Shape of the trapped object. The only demand on 
the trapped object is that it does not rotate in re- 
sponse to forced movement, at least not in a man- 
ner that gives rise to a response in the detection 
system. This condition is fulfilled if the particle 
is spherical, or asymmetric but strongly trapped. 



Due to the height dependence of the drag coeffi- 
cient, even a spherical bead rotates when trans- 
lated close to a surface. However, this rotation dis- 
sipates a negligible amount of energy compared to 
the translational dissipation and hence does 
not affect the calibration. If the trapped object is 
elongated and weakly trapped, it may wobble in 
the trap in response to the oscillating liquid, and 
give rise to detection of false movement by the pho- 
todiode. But this is not a problem with commer- 
cial micro-spheres because they are highly spheri- 
cal, as are many small biological objects, such as 
lipid droplets or vesicles. 

7. Shape of the trapping potential. Throughout this 
paper we assumed a parabolic trapping potential, 
but the method is not limited by this assumption. 
By choosing large drive amplitudes A, it is possible 
to map the shape of the potential and calibrate it. 
This is done by analyzing the higher modes of the 
bead's motion that arise in response to a sinusoidal 
drive in a nonlinear trapping force-field. 

APPENDIX B: HYDRODYNAMICS 

Sometimes it is desirable to work with a flow-cell with 
open ends, e.g. to facilitate the exchange of buffer solu- 
tion. When the ends are open, the liquid between the 
two cover-slips generally does not move with the same 
amplitude or phase as the flow-cell. 

To estimate the size of this effect, we consider the mo- 
tion of a liquid contained between two infinite, paral- 
lel planes which are moved identically, parallel to them- 
selves, in a simple sinusoidal fashion. In an experiment, 
two of the four sides of the flow-cell are sealed, which 
increases the drag force on the liquid. Also, the surface- 
tension at the openings helps to force the liquid to move 
with the flow-cell. Thus, the result below for two infi- 
nite, parallel planes exaggerates the liquid's motion — it 
is a worst-case scenario. 

The velocity of the oscillating planes (the cover-slips) 
is zero in the y and z direction, and 

Vdrivc — ^AcOs{LUt) (Bl) 

in the x-direction, where uj = 27r /drive- The only external 
forces on the liquid are the shear forces arising from the 
no-slip boundary condition between the liquid and the 
accelerated planes. The liquid moves with the acceler- 
ated planes with a time-lag determined by the balance 
between the inertia of the liquid and the shear forces in- 
side the liquid. Thus, we see that the velocity of the 
liquid can be written as u(x, y, z, t) = {u{z, t), 0, 0). 

The equation of motion for the liquid is obtained 
from the linearized Navier-Stokes equations with con- 
stant pressure (see [2^1: 

9u , , d'^u , , 
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where v is the kinematic viscosity of the hquid. The 
relevant solution to this equation can be written 

u{z,t) — a(z) cos(wt) + sin(ti>t) (B3) 

= Aiiquid(^) COs(a;(t - tphasc(2))) , (B4) 



1.000 



where ^uquid = Va^ + b'^ and w^phase = arctan(&/a). 
The coefficients a and b can be found by substituting 
Eq. lUll) into Eq. jB2l) : 



(B5) 



a = ~-b" and 5= -a" 



OJ 



from which we find 



(B6) 



where ' indicates differentiation with respect to z. 
The solution to Eq. ljB6p is a linear combination of 
exp{^^^y^z) whose four coefficients are determined 
by the symmetry requirement a{z) = a{d — z) and the 
boundary conditions a(0) = ijjA and a"(0) = 0: 

z - d/2 z- d/2 
a[z) = ci cos( ) cosh( ) 

,z-d/2^ ,z-d/2, , , 
+ C2sin( -^)sinh( , (B7) 



where 



Cl = LUA 

C2 = toA 



cos(^) cosh(^) 



2C' 



cos2(A)+sinh'(^) 

sin(^)sinh(j^) 
cos2(A)+sinh2(A) ' 



and 



(B8) 
(B9) 

(BIO) 



is the depth of penetration of the shear-wave into the 
fiuid. The expression for b follows from Eq. (jB5|) and is 
identical to Eq. (|B7(I except for ci and C2 swapping places. 
When d^ C the shear-wave's amplitude decreases expo- 
nentially as a function of z. When d ^ C the liquid 
co-moves with the planes as a solid body. 

With the parameters used in the experiment far from 
surfaces, described in Sec.0 i.e., /drive = 28 Hz, the pen- 
etration depth C = 103 fim is larger than the thickness of 
the flow-cell d ~ 60 /im. With these parameters, the am- 
plitude of the liquid's motion differ at most 0.25% from 
the amplitude of the stage's motion, see Fig. 01 This 
result is independent of the amplitude A but is quite 
sensitive to d and /drive- Thus, if working away from the 
surface it is important to choose /drive so that it matches 
the thickness of the cell, i.e., so that d/C < 1. 

It is illustrative to calculate the velocity of the liquid 
midway between the two planes, since this is where the 
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FIG. 4: Hydrodynamic prediction of the response of the liq- 
uid to the motion of the cover-slips. Abscissa: height above 
cover- slip 2. Ordinate: Amplitude of liquid motion yl]iquid(z), 
Eq. 1B4L in units of the amplitude of the oscillating planes 
^stagc = A. Parameters are the same as used in the ex- 
periment far from surfaces (/drive ~ 28 Hz, A — 208 nm, 
r = 23.0°C). 



motion of the liquid differs the most from that of the 
planes: 



m(— , t) = Cl cos{ujt) + C2 sin(a;i) 



(BlI) 



Aiiquid( 2) cos(a;(t - tphase(d/2))) ,(B12) 



where now Auquidi^) = \/ + cj and wtphasc(c?/2) = 
arctan(c2/ci). When the cover-slips are close together 
or the drive frequency is very low ^ ^ 1, and we have 

w(|,i) = Wdrivo(0; i-^-i the fluid co-moves with the plaucs. 
If the cover-slips are far apart or the drive-frequency is 
high Eq. ljB12|l becomes: 



u{-,t) = 2wAe"^ coa{ujt - d/2C), d > 2C , (B13) 

i.e., the amplitude decreases exponentially as d/( grows 
and the motion of the liquid is phase-shifted relative to 
the motion of the planes. These behaviors are illustrated 
in Fig. El From this we see the importance of choosing a 
drive-frequency small enough for the fluid to follow the 
flow-cell in the region where the trap is to be calibrated. 



APPENDIX C: POWER SPECTRAL DENSITY 
FOR FINITE MEASUREMENT TIME 



In Section IIIII we calculated the PSD for the beads' 
motion, assuming inflnite measurement time as well as 
continuous sampling. In a real experiment, data are col- 
lected for a flnite time tmsr and with a flnite sampling fre- 
quency /sample- In what follows we are going to assume 
continuous sampling in time, since the aliasing caused 
by finite sampling time (i) is irrelevant for the discussion 
here, (ii) is easily accounted for 2], and (iii) does not 
occur at all, if over-sampling data acquisition electronics 



10 




FIG. 5: Hydrodynamic prediction of the amplitude of liq- 
uid motion midway between the two cover-slips. Abscissa: 
The distance between two planes d, measured in units of the 
penetration depth Eq. IBIOII . Ordinate: Amplitude of the 
liquid-oscillations midway between the two planes Anquidi^) 
Eq. IIBI2II . in units of the amplitude of the oscillating planes 
^stage = A. Notice the plateau at unity for small values of 
d/(^ where the fluid moves with the planes, and the exponen- 
tial decay to zero at larger values where the liquid no longer 
follows the planes. The values of d/( for the experiments in 
bulk and close to the surface are marked by a square and a 
circle, respectively. 



The recording time ^msr is easily chosen to be an integer 
multiple of the period of the stage's motion /drive- It 
can be done even after recording was finished, simply by 
discarding an incomplete stage period from the recording. 
A simplification is thereby achieved, /drive — fk for some 
integer k, so the spike consists of a single datum 

'^fe = ^^A.,/a.™ (C4) 

and the PSD becomes 

1 j2 ^2 

-Pfe(imsr) — ^2^J2 _|_ J2 _|_ J2 '^/fe. /drive (^5) 

where Sf^j^^^^^ is Kronecker's delta. This is the discrete 
version of the expression given in Eq. 0. We also see 
that the height of the spike depends linearly on the mea- 
surement time tinsr- When the measurement time is in- 
creased the spike approaches Dirac's delta function 

imsr <^/fc, /drive ~* ^if - /drive), tmsr OO (C6) 

and the expression in Eq. ^ is regained. 

APPENDIX D: HYDRODYNAMICALLY 
CORRECT POWER SPECTRUM ANALYSIS 



is used. Readers feeling uneasy about the idea of contin- 
uous sampling need only remind themselves that it is a 
mathematical ideal, a limit that is achieved in practice by 
choosing /Nyq = /sampie/2 > /c and, if ahasing occurs, 
by leaving frequencies with significant aliasing out of the 
analysis. Equations (H5) and (H6) in indicate which 
frequencies have significant aliasing for a given desired 
accuracy. 

For the finite-time continuous-time Fourier transform 
of Eq. (jSj we now have 



Xk 



tmsr/2 



t„r=r/2 



(CI) 



and the expectation value for the one-sided {fk > 0) PSD 
becomes: 



-Pk (^msr ) 



2{\S:k\ ) _ D + ifdrivcrkAy 



where 



rk 



^\Pk+m 



1 Sin(7r(/fc — /drive )^msr) 
i.fk /'drivc)^msr 



(C2) 



(C3) 



gives the shape of the spike in the PSD, referred to as 
leakage 0; see Fig. 5] for peak shapes. Here, fk = 
fc/^msr where k is an integer. We have ignored cross terms 

of the type Sin(7r(/fe - /drive)imsr) Sin(7r(/fc -I- /drive)imsr) 

in Eq. ljC3|) . as they are typically several orders of mag- 
nitude smaller than the terms retained. 



No general statement can be made regarding which ac- 
curacy to expect from a Lorentzian fit. The importance 
of omitted effects depends on several circumstances 
As a rule of thumb, we recommend the use of the hydro- 
dynamically correct theory for Brownian motion, if less 
than 10% systematic error is desired. User-friendly soft- 
ware exists that calibrates with the full theory 0, ^3 j 
and one is left with fewer concerns regarding reliability 
of results. 



1. Hydrodynamically correct power spectrum 

The hydrodynamically correct power-spectrum of 
Brownian motion of a micro-sphere that is trapped by 
a Hookcan force, is 



phydro(^^^/^)^ph.ydro^phyd™ 



response ' 



where 

;jhydro 



P7^™(/, R/£) = 



(Dl) 



(D2) 



(i?o/7r2)Re{7/7o} 



(/c,o + / Im{7/7o} - /V/m,o)2 + (/ Re{7/7o})2 ' 

and 

PrZtonsM,R/e)= (D3) 
KA/drive|7/70|)''^(/- /drive) 

(/c.o + / Ini{7/7o} - /V/rn,o)2 + (/ Re{7/7o})2 ' 
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Equation ljD2|) is the same expression as given in '5, 
Eq. 34] (Note however Here ■y{f,R/£) is the 

frequency-specific drag coefficient that results from solv- 
ing Stokes equations for the micro-sphere undergoing 
linear harmonic motion with frequency /. We assume 
this motion takes place parallel to an infinite planar 
surface that is located a distance £ from the center of 
the microsphere with radius R, and give an approxi- 
mate expression for 7(/, R/£) below. Using the notation 
7o = 7(0, 0) — (STTpvR, where p is the density of the sur- 
rounding liquid and v its kinematic viscosity, the parame- 
ters in Eqs. m2IU3ll are Z^o = fcBT/70, /c,o = K/(27r7o), 
and fm,o — 7o/(27rm) with m the mass of the micro- 
sphere. All three parameters are independent of £. They 
describe bulk values, while all i'-dependence of the PSD 
is found in "f{f, R/£)/^q. In this manner, consistency of 
experimental power spectra that are recorded with the 
same bead and trap at different distances to a surface, 
is very easy to check: The values returned for Dq and 
fmfi should all be the same, while variation in /c_o arise 
from variation due to optical aberrations (that depend 
on £ [HI)- 

2. The drag coefficient 

The drag coefficient 7(/, R/£) is not known exactly for 
R/£ > 0. Deep in bulk, R/£ — 0, and one has Stokes' 
exact result SOil 



7(/,0) 



(D4) 



7stokes(/) = 70 I 1 + (1 - i)W ^ " ^Ij- 



where = vHttR^ ). A t zero frequency one has Faxen's 
approximate result |l3 . IsH , 



7(0,i?/f) = 

7Faxc„(i?/^) 



(D5) 



70 



1 _ iLR. 



45^" 



An exact result exists (32j and has been shown numer- 
ically to agree with Eq. (|D5|) to within 1% for 
£/R>lA. 

An approximate expression for 7(/, R/ £) is given in ^2, 
Eq. (33)]. At zero frequency, it reproduces Faxen's result 
up to and including its first order term in R/£. In the 
spirit of Fade approximants, we rearranged the terms of 



i 



Eq. (33)] to 



(D6) 



7Stokos(/) 



_9__R 

16 i 



1 



1-i 
3 



2i/ 
9/,. 



(1 



-(1-0- 



where 5 — R^J f^j /, and the series in R/ £ in the denomi- 
nator is not known beyond the explicitly shown terms. 
For / = this expression reproduces Faxen's result 



up to and including its second order term in R/ £: For 
£/R > 1.5 it deviates by less than 0.5% from Faxen's re- 
sult, and for £/R > 5 the deviation is smaller than 0.1%. 
This is a significantly better approximation to Faxen's 
result, so it must also be a significant improvement at 
low frequencies, and possibly at all frequencies. 

We found that Eq. (|D2|I agrees with our experimental 
PSDs down to significantly shorter distances from the 
surface when j{f, R/£) in Eq. (|D6|) is used instead of 
Eq. (33) in @], see [Mil. 



3. How to fit tlie ^-dependent power spectrum 

In order to fit the theoretical power spectrum in 
Eqs. ljDlHD3|l to the experimental one, the distance £ 
to the surface must be known. We do know the differ- 
ences between the various distances at which we recorded 
power spectra. It is only the location of the coverslip sur- 
face that is unknown. 

This location can be determined in several ways: From 
the fluorescence induced by an evanescent wave, from the 
effect of the hydrodynamic interaction of the trapped 
bead with the surface, or by analyzing interference or 
diffractive patterns jlj, references therein, and ^2(1 .361]]. 
Yet another method relies on pulling on DNA constructs 
along the axial direction [s^]- Most of these methods 
require a modification of the experimental setup. We 
used the hydrodynamic interaction of the bead with the 
surface, similar to that described in but taking ad- 
vantage of an imp roved understanding of the focal shift, 
as detailed in |ll| . This method requires no extra instru- 
mentation. 

In principle, one might fit Eqs. (|D1IID3|I simultane- 
ously to all power spectra recorded with a given bead 
in the trap, and let the location of the surface be a fit- 
ting parameter. We used a simpler procedure, based on 
the observation that the distance to the surface is needed 
only with moderate precision in order to determine the 
fitting parameters Dq and /c,o with high precision. This 
is because the distance occurs only in Eq. (|D6|I in terms 
that play a subdominant role in Eqs. (|D1IID3|I . 

We fitted Eq. (jD2|l to the thermal background and 
compared the results from inserting either 7 = 7Faxcn 
(turning Eq. (jD2|) into a Lorentzian) or 7 = 7Faxcn x 
7stokcs/7o- The distance-dependent results for /c and D 
were found to agree within 1% for the two approaches. 
This agreement is partly artificial: We described the par- 
asitic filtering in our diode with % Eq. (35)], fitt ing its 
two parameters. This characteristic function absorbed 
most of the difference between 70 and 7stokos(/)- 

Einstein's relation then gives a distance-dependent 
7cx = kBT/D, while Faxen's law should give the dis- 
tance dependence, 7ox ~ 7Faxen {R/£). We therefore fit- 
ted Faxen's law to our ^-dependent values for 7ex, using 
the location of the surface as a fitting parameter. The 
good fits confirmed the procedure OH ■ 

With £ thus known, we know ^{f,R/£) in Eq. (|D6p . 
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and can fit the experimental PSDs with the more cor- 
rect theoretical PSD that results from using Eq. HU6() 
in Eq. I)D2(I . and Dq and fc^ as fitting parameters 
with their original interpretation: Dq ksT/^Q and 
/c,o = K/(27r7o). 

Note that although depends on the radius of the 
bead and on the kinematic viscosity of the surrounding 
fluid, it is systematically correct to insert the manufac- 
turer's value for the radius and the room temperature 
value for the viscosity in f^. This is because in Eq. HU6() 
these small errors occur in a term that is already a small 
correction. 



4. Hydrodynamically correct calibration factor 

The hydrodynamically correct theory described above 
gives 

W^th.hydro = J PXonsedf - (D7) 
if cfl/f drive + Im{7/7o} - /drive//m,o)^ + (Rc{7/7o})2 



This expression does not change the calibration factor /3 
by much, compared to what Eq. would give. This 
is because Einstein's approximate theory for Brownian 
motion that lead to Eq. (|13l) . is a low-frequency approx- 
imation, and /drive is a low frequency compared to the 
characteristic frequency of the frequency dependent 
drag. At room temperature [f dmc/ fuY^"^ — 7.2 -lO""^ for 
the 1.54 /im diameter beads we drove at /drive ~ 28 Hz, 
and ifdi-Wc/fuY^^ = 2.6 • 10"^ for the 528 nm diame- 
ter beads we drove at /drive = 32 Hz. Inertial effects 
are two orders of magnitude smaller 2]. Comparing 
with Eq. (inji, we see that, in bulk, the hydrodynami- 
cally correct theory in our case introduces corrections of 
1.5% and 0.5%, respectively, in the interpretation of the 
power in the spike. In contrast, near a surface the correc- 
tions to Eq. (|13|l are of the opposite sign and a smaller 
magnitude than the corrections in bulk: For £/R = 10, 
/drive = 32 Hz, fcfi ^ 2000 Hz, and 2R = 528 nm the 
correction is —0.25%. 



[1] K. C. Neuman and S. M. Block, Review of Scientific In- [14] 
struments 75, 2787 (2004). 

[2] K. Berg-S0rensen and H. Flyvbjerg, Review of Scientific 

Instruments 75, 594 (2004). [15 

[3] K. Berg-S0rensen et al, Review of Scientific Instruments [16 
77, 063106 (2006). 

[4] S. C. Kuo and M. R Sheetz, Science 260, 232 (1993). 

[5] R. Simmons, J. Finer, S. Chu, and J. Spudich, Biophys. [17 
J. 70, 1813 (1996). 

[6] A common approach calibrates far enough above the cov- 

erslip surface to ensure that uncertainty about the height [18 
does not affect the calibration significantly. This calibra- 
tion is then used closer to the surface. If an oil-immersion [19 
objective is used, this is an error-prone approach: Due to 
the difference in refractive index between oil and water, 
the trap stiffness decreases rapidly with distance from the [20 
surface — 10 fim from the surface, the stiffness is typically 
reduced by more than a factor of two for beads with a [21 
diameter of ~ 500 nm [Tll ligl . I20II . Thus, one detrimental 
effect is swapped for another, and calibration errors can [22 
be large. 

[7] S. B. Smith, Y. Cui, and C. Bustamante, Science 271, [23 
795 (1996). 

[8] M. J. Lang, C. L. Asbury, J. W. Shaevitz, and S. M. [24 

Block, Biophysical Journal 83, 491 (2002). 
[9] K. C. Vermeulen et al., Review of Scientific Instruments [25 
77, 013704 (2006). [26 
[10] M. Capitanio, R. Cicchi, and F. S. Pavone, Eur. Phys. J. 

B. 46, 1 (2005). [27 
[11] E. Schaffer, In preparation (2006). [28 
[12] U. Bockelmann et al, Biophysical Journal 82, 1537 
(2002). 

[13] See Appendix|n|for how to square Dirac's delta-function [29] 
in Eq. @. 



J. Happel and H. Brenner, Low Reynolds number hydro- 
dynamics (Noordhoff International Publishing, Leyden, 
1973). 

F. J. Harris, Proc. IEEE 66, 51 (1978). 
I. M. Tolic-N0rrelykke, K. Berg-S0rensen, and H. Flyv- 
bjerg, Computer Physics Communications 159, 225 
(2004). 

P. M. Hansen, I. M. Tolic-N0rrelykke, H. Flyvbjerg, and 
K. Berg-S0rensen, Computer Physics Communications 
174, 518 (2006). 

S. Hell, G. Reiner, C. Cremer, and E. H. K. Stelzer, Jour- 
nal of Microscopy-Oxford 169, 391 (1993). 
E.-L. Florin, A. Pralle, E. H. K. Stelzer, and J. K. H. 
Horber, Applied Physics A: Materials Science & Process- 
ing 66, 75 (1998). 

K. C. Neuman, E. A. Abbondanzieri, and S. M. Block, 
Optics letters 30, 1318 (2005). 

E. J. G. Peterman, F. Gittes, and C. F. Schmidt, Bio- 
physical Journal 84, 1308 (2003). 

W. Denk and W. W. Webb, Applied Optics 29, 2382 
(1990). 

K. Svoboda, C. F. Schmidt, B. J. Schnapp, and S. M. 
Block, Nature 365, 721 (1993). 

A. Pralle et al., Microscopy Research and Technique 44, 
378 (1999). 

J. Howard and A. J. Hudspeth, Neuron 1, 189 (1988). 
K. Berg-S0rensen, L. Oddershede, E.-L. Florin, and H. 
Flyvbjerg, Journal of Applied Physics 93, 3167 (2003). 
A. Rohrbach, Physical Review Letters 95, 168102 (2005). 
L. D. Landau and E. M. Lifshitz, Fluid Mechanics, Vol. 6 
of Course of Theoretical Physics (Pergamon Press, Ox- 
ford, 1959). 

Equation (34) in 2] contains a typographical error: The 
parameter fm in this equation is fm,o = 7o/(27rm), 



13 



i.e., not the same as the parameter fm, which occurs 
in Eq. (32) in |3. Systematic naming is achieved, if 
the latter is referred to as fm*,o ~ 7o/(27rm*), i.e., 
only its name is changed, not its definition in terms of 
m* — m -\- /Z, where p is the mass density of the 

liquid |39| . This error is not repeated in the programs 
described in H^ITl. 

[30] G. G. Stokes, Trans. Cambridge Philos. Soc. 9, 8 (1851). 

[31] H. Faxen, Ph.D. thesis, Uppsala, 1921. 

[32] M. O'Neill, Mathematika 11, 67 (1964). 

[33] A. J. Goldman, R. G. Cox, and H. Brenner, Chem. Eng. 
Sci. 22, 637 (1967). 

[34] Consequently, Eq. <D6ll 's expression for 7(f , Rll) has 
been introduced in the programs described in Ua . [7^ 
in the latest versions found in Computer Physics Com- 
munication's program library |40|. 



[35] Note, however, that a function can approximate another 
function well, even when their derivatives do not not ap- 
proximate each other well. This is the case here: If one 
considers how Faxen's result "/{0,R/£) is approached by 
^{f,R/t) by expanding Eq. IID6t in powers of \J j I fv^ 
the approach is incorrect to first order. 

[36] P. M. Hansen, J. K. Dreyer, J. Ferkinghoff-Borg, and 
L. Oddershede, Journal of Colloid and Interface Science 
287, 561 (2005). 

[37] C. Deufel and M. D. Wang, Biophys. J. 90, 657 (2006). 

[38] M. D. Wang et ai, Biophys. J. 72, 1335 (1997). 

[39] K. Berg-S0rensen and H. Flyvbjerg, Erratum to 0] (To 
appear) . 

[40] K. Berg-S0rensen and H. Flyvbjerg, private communica- 
tion. 



